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Abstract. We present the first near-infrared A"-band long-baseline interferometric measurement of the nucleus of the prototype 
Seyfert2 Galaxy NGC 1068 with resolution A/B ~ 10 mas obtained with the Very Large Telescope Interferometer (VLTI) and 
the two 8.2 m diameter Unit Telescopes UT2 and UT3. The adaptive optics system MACAO (Multi Application Curvature 
Adaptive Optics) was employed to deliver wavefront-corrected beams to the A"-band commissioning instrument VINCI. A 
squared visibility amplitude of 16.3 ±4.3% was measured for NGC 1068 at a sky-projected baseline length of 45.8 m and 
azimuth angle 44.9 deg. This value corresponds to a FWHM of the A"-band intensity distribution of 5.0 ± 0.5 mas (0.4 ± 0.04 pc 
at the distance of NGC 1068) if it consists of a single Gaussian component. Taking into account A"-band speckle interferometry 
observations (Wittkowski et al. 1998 . Weinberger et al. 1999 . Weigelt et al. 2004 1, we favor a multi-component model for the 
intensity distribution where a part of the flux originates from scales clearly smaller than ~ 5 mas (< 0.4 pc), and another part 
of the flux from larger scales. The A"-band emission from the small (< 5 mas) scales might arise from substructure of the dusty 
nuclear torus, or directly from the central accretion flow viewed through only moderate extinction. 

Key words. Techniques: interferometric - galaxies: nuclei - galaxies: Seyfert - galaxies: individual: NGC 1068 



1. Introduction 

The Seyfert galaxy NGC 1068 harbors one of the brightest and 
closest active galactic nuclei (AGN). Active galaxies appear as 
types 1 and 2, where the spectra of the former exhibit broad and 
narrow emission lines, and those of the latter show only narrow 
lines. Antonucci & Miller ( 1985 1 suggested that the broad-line 
emission regions are located inside an optically and geometri- 
cally thick disk and that central continuum and broad-line pho- 
tons are scattered into the line-of-sight by free electrons above 
and below the disk. Depending on the observer's viewing an- 
gle, the broad-line emission region is either obscured or not. 
This suggestion is now widely accepted and has evolved to the 
so-called "unified scheme of AGN" (e.g., Antonucci] 1993 1. 

Various theoretical models of the postulated dusty tori have 
been presented by, for instance, Krolik & Begelmann (1988 1, 
Pier & Krolik (11^921 . Granato & Danese (JT994}, Efstafhiou 
& Rowan Robinson (1995 1, Manske et al. (19981, Nenkova 
et al. d2"0"0"2l . and Vollmer et al. J2004> . Owing to the lack of 
spatially resolved observations, the models were usually com- 
pared to the integrated spectrum of the core of NGC 1068. It 
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turned out that torus models with a wide variety of geometries, 
spatial extensions, and optical depths are consistent with this 
spectrum. In addition, the nature of the central emission source 
itself could not be well constrained because the flux spectrum 
of the very inner nuclear engine could not be separated from 
the emission of surrounding material. 

Several high-resolution infrared observations of NGC 1068 
showing a compact central IR core and surrounding structure 
were carried out by, for instance, Thatte et al. ([1997 ), Rouan et 
al. i fT9l?8l Wittkowski et al. i fT9l?8l Weinberger et al. ifTW9l, 
and Bock et al. .2000. . Wittkowski et al. {E>98 1 presented both 
a AT-band visibility function up to spatial frequencies corre- 
sponding to a baseline of 6 m, as well as the first 76 mas res- 
olution AT-band image of the nucleus of NGC 1068 obtained 
by bispectrum speckle interferometry. The compact central IR 
core was resolved with a Gaussian FWHM of ~ 30 mas (2pc). 
New AT-band and //-band bispectrum speckle interferometry 
observations were performed by Weigelt et al. (2004). 

Very recently, the first interferometers consisting of 8-10 m 
class telescopes started operations, and they have already suc- 
ceeded in observing AGN with much higher spatial resolu- 
tions. The potential of optical/infrared interferometry to inves- 



tigate AGN was recently discussed by Wittkowski et al. ( 2003 1. 
Swain et al. ( 2003 1 reported the first AT-band interferometric ob- 
servations of the Seyfert 1 galaxy NGC 4151 obtained with the 
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Table 1. Calibration sequence of NGC 1068. JD is the Julian date of the observation, N the number of processed interferograms, B 
and Az the projected baseline length and azimuth angle (east of north), and /j 2 the obtained coherence factors with their statistical 
errors. The interferometric efficiency (IE) given in italic characters is the value adopted for the calibration of the NGC 1068 data. 



JD 


N 


B(m) 


Az 


/i 2 + stat. 


IE ± stat. ± syst. 


V 2 ± stat. ± syst. 




Target 


2452948.717 


250 


45.79 


44.87 


0.054 + 0.017 


0.340 ± 0.009 ± 


0.001 


0.158 +0.049 ± 


0.001 


NGC 1068 


2452948.722 


214 


46.00 


45.13 


0.062 + 0.031 


0.340 + 0.009 ± 


0.001 


0.182 ±0.091 ± 


0.001 


NGC 1068 


2452948.761 


407 


46.62 


45.40 


0.319 + 0.005 


0.348 + 0.006 + 


0.001 






HD 20356 


2452948.767 


438 


46.64 


45.62 


0.304 ± 0.005 


0.331 ±0.005 + 


0.001 






HD 20356 


2452948.772 


416 


46.63 


45.80 


0.315+0.006 


0.344 + 0.007 ± 


0.001 






HD 20356 



85 m baseline of the Keck Interferometer. Jaffe et al. ( 2003) re- 
ported on the first mid-infrared interferometric observation of 
the dusty torus of NGC 1068 using VLTI/MIDI. 

In the present letter, we report on the first /f-band long- 
baseline interferometric observation of NGC 1068. 



2. Observations and data reduction 

The NGC 1068 interferometric data were obtained with the 
Very Large Telescope Interferometer (VLTI) and the K- 
band commissioning instrument VINCI (Kervella et al. 2000 
Kervella et al. 2003 ), used with the fiber-based beam combiner 
MONA, on Nov. 4, 2003. The UT2-UT3 baseline with 47 m 
ground length was used. Both telescopes were equipped with 
the Multi Application Curvature Adaptive Optics (MACAO) 
system. These data were taken in the framework of the com- 
missioning of MACAO- VLTI. 

MACAO systems. For a detailed description of the MACAO 
systems see Arsenault et al. (2003). The MACAO systems are 
four curvature adaptive optics systems installed at the Coude 
focus of each UT (the systems for UTs 1 and 4 will be installed 
in the fall of 2004 and first quarter of 2005). The Coude train 
mirror M 8 is replaced by a 60 actuator bimorph deformable 
mirror in a tip-tilt mount. The mirror M 9 is a dichroic which 
transmits visible light to the wavefront sensor and reflects the 
wavelength range from 1 to 13//m to the VLTI recombina- 
tion laboratory. During commissioning, the two MACAO sys- 
tems on UTs 2 and 3 showed very similar and consistent re- 
sults. Strehl ratios of 65% were obtained on bright guide stars 
(V < 11) under optical seeing conditions of up to 0.8". On 
fainter sources with V magnitudes of 14.5 and 16, Strehl ra- 
tios of 45% (seeing 0.60") and 25% (seeing 0.55") were ob- 
tained, respectively. For use with the single-mode fiber instru- 
ment VINCI, a high Strehl ratio ensures a high concentration 
of light onto the 56 mas diameter fiber core of VINCI, hence 
allowing the interferometric observation of faint sources such 
as NGC 1068. For both MACAO- VLTI systems the reference 
source was the nucleus of NGC 1068 itself, and the same con- 
trol parameters were used for both AO systems. No neutral den- 
sity filters were used, the main loop gain was 0.50. The optical 
seeing at the time of observation was ~0.9". It is difficult to 
estimate the MACAO performance for our NGC 1068 obser- 
vations since the source is extended and an acquisition image 
could not be taken. 
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Fig. 1. Average wavelets power spectral density (WPSD) of the 
464 processed interferograms of NGC 1068. The subtraction of 
the background noise from the processed fringes' power peak 
left no residual bias on the final WPSD. The power integration 
is done between wave numbers 2000 and 8000 cm -1 . 



Interferometric data. 1 000 interferograms of the nucleus of 
NGC 1068 were obtained in two series of 500 scans. The 
fringe frequency was 216 Hz, a compromise between faint ob- 
ject sensitivity and immunity to the atmospheric piston effect. 
Of the recorded interferograms, 464 were processed success- 
fully by the VINCI data reduction software as described by 
Kervella et al. (2004). Table [2 lists the observational details 
for the NGC 1068 and the calibrator data, as well as the re- 
sulting visibility values. Figure ^ shows the background- and 
noise-corrected average wavelet power spectral density of the 
NGC 1068 interferograms. The fringes' power peak around 
wave number 4500 cm -1 is not affected by any significant 
power spectrum bias, despite the faintness of the source. The 
effective wavelength of the NGC 1068 observations is approx- 
imately 2.18jum. The K5 giant HD 20356 (0 VD = 1.81 + 
0.02 mas) from Cohen et al. (fT999l and Borde et al. d2002l was 
used as calibration star (effective wavelength A = 2.181 fim). 
Observational parameters and data reduction of the calibration 
star were identical to those of NGC 1068. The systematic error 
induced by the calibrator on the final visibility values is neg- 
ligible compared to the statistical error. A calibrated squared 
visibility of V 2 = 16.3 + 4.3 % for NGC 1068 was obtained for 
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a sky-projected baseline with length B = 45.839 m and azimuth 
angle Az = 44.93 deg (east of north). 

Photometric estimate. Any photometric estimates using a 
single-mode fiber instrument is, in general, difficult due to the 
large and rapid fluctuations of the coupling of the object light 
into the fiber core. However, the MACAO systems keep a large 
fraction of the object light inside the Airy disk, stabilize the 
injected flux, and a photometric estimate can be attempted. 
For commissioning purposes, a number of stars of various K- 
band magnitudes were observed on Nov. 3-5, 2003 with the 
same VLTI configuration as used for our NGC 1068 observa- 
tions. The relation between the observed flux values and the 
/f-band magnitudes is consistent with the expected expo- 
nential function, so that the attempt of an absolute photomet- 
ric calibration is reasonable. The best-fit relations between mg- 
and the photometric fluxes Pa and Pg (in ADU/s) are Pa = 
3.50 10 6 exp(-0.906 m K ) and P B = 1.46 10 6 exp(-0.901 m K ). 
The residual dispersions on are <ta = 0.5 mag and cr B = 
0.6 mag. 

Since the performance, i.e. the Strehl ratio, of the MACAO 
systems for our extended source NGC 1068 is unknown, but 
very likely lower than that of the single bright stars (see above), 
the application of this photometric estimate on NGC 1068 can 
only give us a lower limit of the NGC 1068 flux in our field of 
view (FOV) of F K > 130 ± 60mJy (m K < 9.2 + 0.4). The size 
of the K-bmd Airy disk of the UTs with MACAO is 56 mas. 
In the absence of atmospheric turbulence, the FOV of our mea- 
surements would correspond exactly to the projection of the 
single-mode fiber on the sky. As this mode is matched by de- 
sign to the diffraction pattern of the UTs, the effective FWHM 
of our FOV would also be 56 mas. In practice, however our 
FOV is slightly larger due to the presence of residual uncor- 
rected speckles, i.e. the limited Strehl ratio. 



3. Discussion 

Intensity distribution of the /f-band emission. Our mea- 
sured VLTl/VINCI ^T-band squared visibility amplitude of 
\V\ 2 = 16.3 + 4.3 % at a projected baseline length of B = 
45.8 m corresponds to a single-component Gaussian inten- 
sity distribution with FWHM 5.0 ± 0.5 mas (taking into ac- 
count the VINCI broad-band K filter). With the distance to 
NGC 1068 of 14.4Mpc (Bland-Hawthorn et al. IT9971 . this 
corresponds to a linear FWHM of 0.4 + 0.04 pc. This means 
that the VLTl/VINCI visibility value is consistent with a multi- 
component model only if at least 40% (|V| = 40.4 %) of the 
flux (i.e. > 50mJy) originates from scales <5mas (Gaussian 
FWHM < 6.7 mas or < 0.5 1 pc at the 3 cr level). 

The K-bwd visibility values obtained by speckle inter- 
ferometry at spatial frequencies up to a baseline of 6-10 m 
(Wittkowski et al. 1998 SAO 6 m telescope; Weinberger et 
al. [T9991 Keck 10 m; Weigelt et al. 120041 SAO 6 m) are 
consistent with a single-component Gaussian intensity dis- 
tribution with azimuthally averaged FWHM ~ 30 mas (up to 
B = 6m; Wittkowski et al.[1998), and also with a larger ~30- 
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Fig. 2. Synthetic squared visibility function for (dashed line) 
a 5 mas Gaussian, (dashed-dotted line) a 30 mas Gaussian, and 
(solid lines) two examples for a two-component model where 
(upper line/ lower line) 56%/ 40% of the total flux comes from 
a 3 mas/ 0. 1 mas Gaussian and the remaining 44%/ 60% from a 
54 mas/ 42 mas Gaussian (FWHM values). Measurements are 
available for baselines B — 46 m reported in this letter (marked 
by the x-symbol with error bar) and up to B = 6 - 10 m 
(Wittkowski et al. 119981 Weinberger et al. 1999 Weigelt et al. 
2004 1. Both measurements are only matched if the small com- 
ponent has a size clearly below ~ 5 mas and if only part of the 
total flux in our FOV arises from this < 5 mas component, and 
another part from a larger component. It is assumed that the 
FOV is the same for all shown models. 



50 mas component plus a much smaller (unresolved) compo- 
nent (Weinberger et al. 119991 Weigelt et al. l2004t . 

The comparison of the visibility measurements obtained by 
these different methods is difficult since the visibility scales 
with the total observed flux in each FOV. However, since the 
speckle measurements show a structure with FWHM -SO- 
SO mas, and the VLTl/VINCI FOV is -56 mas, the total flux 
observed by VLTl/VINCI is very similar to that of the compact 
~ 30-50 mas speckle component. 

Both, the VLTl/VINCI and the speckle measurements 
are consistent with a multi-component intensity distribution 
where > 50 mJy originate from scales < 5 mas or < 0.4 pc 
(VLTl/VINCI), and another part of the flux from larger scales 
of the order of 40 mas or 3pc (speckle). Fig.[2]shows the visi- 
bility models for a 5 mas Gaussian matching our VLTl/VINCI 
measurement, a 30 mas Gaussian matching the speckle mea- 
surements up to B = 6m, and two examples for two-component 
models matching both, VLTl/VINCI and speckle measure- 
ments. It illustrates that two-component models are consistent 
with both measurements only if the small component has a size 
clearly below ~ 5 mas and a flux contribution of clearly less 
than the total flux in our FOV. 

In the following, we discuss the possible origin of the newly 
constrained very compact (< 5 mas) K-bwd component. 

Sources of very compact /f-band emission. Dust at the inner 
cavity of the torus is heated to near the sublimation tempera- 
ture of about 1000-1500K and emits thermal /^-band photons. 
These photons could escape in a direction other than the equa- 
torial plane into our line-of-sight. Depending on the adopted 
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dust properties and luminosity of the central source, the subli- 
mation radius could be as small as about 0.3 pc (Barvainis et al. 
1987 Thatte et al. 1997 1, or may as well be of the order of 1 pc 
or above (Dopita et al. 1998 1. Hence, the size of the inner dust 
cavity seems to be larger than our observed scale of FWHM 
clearly below ~0.4pc. It may more likely coincide with the 
scale of the ~40 mas (~3pc) speckle component. However, 
since the work of Krolik & Begelmann ( 1988 1 it has been dis- 
cussed that the dusty torus very likely does not have a smooth 
uniform shape but may consist of a large number of clumps 
(Nenkovaatal. 120021 Vollmeret al. EOM . Vollmer et al. (2004i 
give a theoretical estimate for the radius of such clumps of 
rci ~ 0.1 pc. Since this size is consistent with our observed 
scale (<0.4pc), one might speculate that light from distinct 
clumps may contribute to our observed flux. Such dust clumps 
can also exist in polar direction at similar distances from the 
nucleus. Another possibility is that free electrons located above 
and below the broad-line emission region in the ionization cone 
at distances of only a fraction of a parsec scatter central light 
into our line-of-sight. 

It has already been discussed by Wittkowski et al. (1998 1 
that a fraction of the /f-band photons originating from the cen- 
tral engine could reach us directly through only moderate ex- 
tinction in the near-infrared, despite the Seyfert 2 type of this 
AGN. The A v could be as low as ~ 10 m (Bailey et al. 1988 1, 
corresponding to Ak ~1.2m assuming standard galactic ex- 
tinction. With the concept of a clumpy torus, the chance of low 
extinction towards the central source is even larger than for a 
smooth dust distribution. In this case, an analysis of the sepa- 
rated flux spectrum of only the very compact (<5mas) com- 
ponent should reveal a type 1 spectrum, which is supported 
by the possible detection of a broad Bry line by Gratadour et 
al. l|2p03 1. Non-negligible /f-band flux contributions of the or- 
der of a few hundred milli-Jansky could arise from the central 
accretion flow (cf. Wittkowski et al. 1998, Beckert & Duschl 
2002 Weigelt et al. 2004 1. Stars are very unlikely to contribute 
significantly to the A'-band flux on scales <0.4pc (see e.g., 
Thatte et al. lT997l 

Comparison to VLTI/MIDI NGC 1068 and Keck NGC 4151 
observations. It is striking that the two NGC 1068 Zf-band 
scales found by our VLTI/VINCI observations and by speckle 
interferometry (Wittkowski et al. 119981 Weigelt et al. l2004> of 
<0.4pc and about 2 x 4pc, respectively, are very similar to 
those of the two NGC 1068 dust components reported by Jaffe 
et al. J2004> based on the V-band VLTI/MIDI observations. 
The latter include a hot (T > 800 K) 0.8 x (< 1) pc and a warm 
(T ~ 320 K) ~ 2.5 x 4pc component. Hot thermal emission 
from inner substructure of the dusty torus, as well as direct light 
from the central accretion flow could explain the ^-band as 
well as the V-band very compact sub-parsec component. It is 
unlikely, though, that the warm VLTI/MIDI component and the 
2 X 4 pc /f-band component can be explained by the same dust 
component. 

The Keck Interferometer observations of the Seyfert 1 nu- 
cleus of NGC 4151 (Swain et al. 2003 1 showed that the major- 
ity of the central (~ 3 pc) /T-band light arises from very small 



scales of < 0.1 pc, probably from the central accretion disk. Our 
/f-band observations of the Seyfert 2 galaxy NGC 1068 show 
that only part of the central light arises from compact < 0.4 pc 
scales while another part arises from larger scales. This differ- 
ence is in line with the unified scheme which predicts that the 
central engine of Seyfert 2 cores is obscured, which leads to a 
larger (up to 100%) relative flux contribution from surrounding 
material. 

4. Conclusion 

We have obtained /f-band interferometric observations of the 
nucleus of the Seyfert 2 galaxy NGC 1068, and hereby show 
that the correlated magnitude of this AGN up to a baseline of 
at least ~ 50 m is such that it can be studied with the VLTI at 
near-infrared wavelengths. 

Together with other observations, we conclude that a K- 
band flux of >50mJy originates from scales clearly smaller 
than about 5 mas or 0.4 pc and another part of the flux from 
larger scales. Our VLTI/VINCI measurement alone sets an up- 
per limit of < 6.7 mas at the 3 <x level to the Gaussian FWHM 
of the very compact component. The origin of this newly con- 
strained small-scale emission can be interpreted as substructure 
of the dusty torus, as for instance part of a clumpy inner cav- 
ity or distinct clumps forming the torus, or as direct emission 
from the central accretion flow viewed through only moderate 
extinction in the K-bwd. 

Acknowledgements. The VLTI observations reported here were made 
possible through the efforts of the whole ESO VLTI and MACAO 
commissioning teams. 
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